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The angular distribution of low-frequency radiation after single scattering 
by an isotropic distribution of relativistic electrons considerably differs from 
the Rayleigh angular function. In particular, the scattering by an ensemble 
of ultra-relativistic electrons obeys the law p — 1 — cos a, where a is the 
scattering angle; hence photons are preferentially scattered backwards. We 
discuss some consequences of this fact for astrophysical problems. We show 
that a hot electron-scattering atmosphere is more reflective than a cold one: 
the fraction of incident photons which become reflected having suffered a 
single scattering event can be larger by up to 50 per cent in the former case. 
This should affect the photon exchange between cold accretion disks and 
hot coronae or ADAF flows in the vicinity of relativistic compact objects, as 
well as the rate of cooling (through multiple inverse- Compton scattering of 
seed photons supphed from outside) of optically thick clouds of relativistic 
electrons in compact radiosources. The backward scattering asymmetry also 
causes spatial diffusion of photons to proceed slower in hot plasma than in 
cold one, which is important for the shapes of Comptonization spectra and 
time delays in the detection of soft and hard radiation from variable X-ray 
sources. 
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1 INTRODUCTION 



In our recent paper (Sazonov, Sunyaev, 1999a) we have called attention to 
the fact that the differential cross-section for Thomson scattering of low- 
frequency photons averaged over an isotropic distribution of relativistic elec- 
trons is substantially different from the Rayleigh phase function, which cor- 
responds to the scattering by cold electrons. The resulting angular function 
is backward-oriented, i.e. photons tend to be scattered backwards, rather 
than forwards (see Fig. 1). It should be immediately emphasized that this 
angular function is a characteristic of the scattering by an ensemble of elec- 
trons. Similar terminology is used when treating the Compton scattering in 
dense plasmas, where collective effects are important (see, e.g., Bekefi, 1966). 

The phenomenon mentioned above results from the combined operation 
of two effects. One is that a photon is more likely to suffer a scattering 
from an electron that is moving toward it, rather than away from it (the 
probability of scattering is proportional to the factor 1 — cos9v/c, where 
V is the electron velocity and 9 is the angle through which the photon and 
electron encounter) . The other effect is that photons emerge after scattering 
preferentially in the direction of the motion of the relativistic electron. The 
angular function in the considered case contrasts with the forward-oriented 
Klein-Nishina angular function which corresponds to the case of scattering 
of energetic photons on a resting electron. 
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In the present paper, we consider some astrophysical implications of this 
pecuhar scattering behaviour of hot plasma. The most evident consequence 
of the backward scattering law is that a larger (than in the case of cold 
plasma) fraction of low-frequency photons incident on the surface of an op- 
tically thick cloud of relativistic electrons will be reflected from it after a 
single scattering by an electron. Therefore, less photons will participate in 
the Comptonization process (through multiple scatterings) inside the cloud, 
and smaller will be the rate at which the electrons lose their energy. We 
derive (in Section 3) , using classical results of the theory of radiative transfer 
in scattering atmospheres, simple formulae for the albedos of atmospheres 
consisting of either mildly or ultra relativistic thermal electrons. The results 
obtained can be useful, in particular, for calculations of the photon exchange 
between cold accretion disks and hot coronae (or outflows) above them, or 
in the study of the illumination of an ADAF accretion flow by external low- 
frequency radiation. 

Another consequence of the backward scattering asymmetry is that the 
coefficient of spatial diffusion of photons in hot plasma should be different 
from that in the case of non-relativistic plasma. This affects the shapes of 
Comptonization spectra and time delays between the detections of soft and 
hard radiation from variable X-ray sources. We discuss this issue in Section 4. 

It is necessary to mention that the phenomenon of anisotropic (backward) 
scattering of low-frequency radiation in a hot plasma has been previously 
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discussed in literature, and its various astrophysical consequences have been 
studied in detail, mainly in the context of the mechanisms of formation of 
hard X-ray spectra in compact sources (Ghisellini et al., 1991; Titarchuk 
1994; Stern et al., 1995; Poutanen, Svensson, 1995; Gierlinski et al., 1997; 
Gierlinski et al., 1999). In particular, the scattering angular function for an 
ensemble of Maxwellian electrons was investigated in the paper by Haardt 
(1993), where approximate formulae were obtained for it based on the results 
of Monte-Carlo simulations. The main purpose of the present paper is to 
present a set of exact analytical formulae for the scattering angular function 
and some related quantities, which are applicable in the ultra-relativistic and 
mildly relativistic limits. 

2 SCATTERING ANGULAR FUNCTION 

The scattering angular function p{fis) gives the probability of scattering 
through a given angle a = arccos/ig of a photon by an ensemble of elec- 
trons. Below p{fis) is normalized so that the mean photon free path in the 
plasma is 

where is the electron number density and ctt is the Thomson scattering 
cross-section. Note that according to the definition J p(yUs) dfis/2 < 1 in 
the general case (because of Klein- Nishina corrections). 

We shall restrict ourselves in this paper to the case in which low-frequency 

5 



radiation is scattered by electrons that obey a relativistic Maxwellian dis- 
tribution. Two opposite limits within this case can be considered: (A) of 
ultra-relativistic electrons, kT^ ^ meC^, and (B) of weakly relativistic elec- 
trons, kTc <^ meC^ (where is the electron temperature). In case (A), the 
following asymptotic formula is applicable (Sazonov, Sunyaev, 1999a): 

p(/is) = 1 - /is- (2) 

This expression exactly corresponds to the ideal case where Te — oo and 
u ^ {u being the initial photon energy). The angular function is backward- 
oriented in this case (see Fig. 1). It is possible to add first-order correction 
terms, one allowing for Klein-Nishina corrections and another taking into 
account the finite energy of electrons, to equation (0) (Sazonov, Sunyaev, 
1999a): 

,^ , ^ hiy kT, . ,2 -1 +3 ln4 + 6r(0, mecV A;re) /mec2 

p(Te, Z/,/is) = l-yUs-6 ^ -{l-Hs) + 



rrieC^ nieC^ 8 \ kTe y 

(3) 

where r{a,z) = x"~^e~^dx is the incomplete Gamma function. Formula 
d^) is a good approximation if kTg ^ 2meC^, hvkT^ ^ 0.01 (meC^)^. 

In case (B), the angular function can be described by the formula (Sazonov 
& Sunyaev 1999b) 

3 r hv kT 
p(Te, z/, /is) = - 1 + ^^l - 2(1 - /is)(l + ^ + 2(1 - 2/is - + 2^^l)^ 

+ (1 - /is)'(4 + 3/^2) ^ + (1 - /i3)(-7 + 14/Xs + 9/^2 - 10/xf) \ 



( 



) 
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+ (-7 + + 9/i,2 - 38^^ + 20/if) 



+ ... 



which is a good approximation when kT^ <^ O.OSmeC^ and hu <^ O.lmgC^. 
The main term in the power series above, 3(1 + /is)/4, is the usual Rayleigh 
function, which corresponds to the non-relativistic case. 

Various examples of the angular function corresponding to the scattering 
on high-temperature electrons, as resulted from Monte-Carlo simulations or 
calculated from equations @ and (^), are presented in Fig. 1. 



3 REFLECTION FROM A HOT PLANE- 
PARALLEL ATMOSPHERE 



The angular laws of scattering by ensembles of Maxwellian electrons de- 
scribed above are readily applicable to any classical problem of radiation 
transfer in scattering atmospheres. In fact, as long as the evolution of the 
photon energy in time and space is not important, i.e. the scattering angular 
function can be considered unchanging (which takes place, e.g., if hu remains 
so small that Klein-Nishina corrections can be ignored), it is possible to ap- 
ply some classical results (found, e.g., in the textbooks by Chandrasekhar, 
1950 and Sobolev, 1963) to hot electron-scattering atmospheres just by sub- 
stituting the ensemble-averaged angular function for the phase function in 
the corresponding formulae. We shall consider below the problem about the 
reflection of light from such an atmosphere. 



Since Compton scattering involves the transfer of energy between elec- 
trons and radiation, it is convenient to carry out the treatment in terms of 
number of photons, rather than in terms of intensity. We aim to determine 
what fraction of photons incident on the atmosphere become reflected having 
suffered only one scattering event. At the same time, that will immediately 
tell us what fraction of the incident photons are capable of penetrating inside 
the atmosphere, undergoing multiple scatterings and thus taking part in the 
Comptonization process (here we are talking about an atmosphere with an 
optical depth r ;^ 1, otherwise a large fraction of photons will traverse the 
atmosphere un-scattered) . Throughout our discussion below, we restrict our- 
selves to the case in which the incident photons are of sufficiently low energy 
so that Klein-Nishina corrections are not important. It should be specially 
noted that in this (Thomson) limit / p{iis) d/ig/'^ — 1- 

In the case of a semi-infinite atmosphere of optical depth r — > oo, the 
number of photons which have suffered a single scattering process, emergent 
in a specified direction n = cos 6, (p (relative to the normal to the atmo- 
sphere), is given by the relation (Chandrasekhar, 1950) 

where /xq, <^o define the direction from which the radiation is incident on the 
atmosphere, Fq is the incident photon flux per unit surface area, and p is the 
scattering phase function. 

We may introduce the notion of atmosphere albedo with respect to singly 
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scattered photons: 

A^^\^^) - / / ^^^HAt) 0; /^O) 0o)/^ 

which is a function of the incidence angle fiQ. 

The calculation in equation (^) is straightforward for two scattering laws 
which are of interest to us: (A) p = 1 — fig and (B) p = 3(1 + /is)/4 (the 
Rayleigh law), for which the single-scattering albedos are, respectively. 



(A)A«(/io) = I 



Z /io 



[B)A«(/^o) = ^ 



- + -/io - /io - -/^o + 3/io - /^o(3 - 2/io + 3/io) In 

O Z Z /io 



Case (B) corresponds to a situation in which the electrons are cold [kT^ <^ 
meC^), while case (A) will be realized, as follows from the results of the 
previous section, if the plasma is very hot {kTe ^ meC^)- The reflection 
properties of the atmosphere that obeys the scattering law p = 1 — /is are 
noticeably different from those of the cold atmosphere, as demonstrated by 
Fig. 2(a). In particular, the single-scattering albedo for normally incident 
photons {fiQ = 1) is 0.25 for the atmosphere consisting of ultra-relativistic 
electrons, whereas the corresponding value for case (B) is 0.168. This means 
that 50 per cent more photons are reflected, having experienced a single 
scattering event, from the hot atmosphere than from the cold one. It is 
evident from Fig. 2(a) that for kT^ ^ 2mcC^ the albedo approaches the value 
corresponding to the limiting, ultra-relativistic case (A). 



If the electron temperature is only mildly relativistic, the appropriate 
expression for the albedo can be found by integrating the angular function 



3 r 8 1 3 1 I 
^^'^(-"o) = o + 9-^0 - -"0 - o-^o + 3/io - /^o(3 - 2iil + In ^ + [-2/io + 



+ - ^/i^ - + 20/i[; - /io(-2 + - 42/i^ + 20/i^) In 

d /io 

where we have quoted only the correction term proportional to kT^/rrieC^. 
This formula is a good approximation if kT^, ^ 0.2mcC^. As follows from both 
equation and Fig. 2a, the single-scattering albedo of the atmosphere con- 
sisting of mildly relativistic electrons is only slightly larger than that of the 
atmosphere with the Rayleigh scattering law. The relative difference is less 
than 5 per cent if kT^ <^ O.lmgC^. We may conclude that in most practical 
applications, it should be sufficiently accurate to use the albedo correspond- 
ing to the cold atmosphere when reflection of radiation from plasmas with 
kTe ^ OArricC^ is investigated. 

Although formulae (|) and are, strictly, applicable only to an atmo- 
sphere whose optical thickness is infinite, a similar increase in the albedo dur- 
ing the transition from the cold to hot electrons occurs in the case where the 
atmosphere is transparent. As follows from results of Monte-Carlo simula- 
tions, the single-scattering albedo with respect to normally incident photons 
is increased by approximately 50 per cent on going from the cold to ultra- 
relativistic electrons for arbitrary values of the atmosphere optical depth. 
An example in which r = 0.5, a value typical for the advection fiows near 

10 



kT, 
rripC^ 



(9) 



accreting compact objects, is presented in Fig. 2(b). 

We can also calculate the single-scattering albedo of an optically thick 
atmosphere for the case of isotropic incident radiation (when /o(a*o, 0o) = 

const): 

(A) >1« = « 0.282, (10) 

(B) 4i = H^Z^ « 0.208. (U) 

We obtain the result that 36 per cent more photons are reflected (after a 
single scattering) from the ultra-hot atmosphere than from the cold one. 

4 PHOTON DISTRIBUTION OVER THE 
ESCAPE TIME FROM A CLOUD OF 
ELECTRONS 

In the previous section we showed that the fraction of low-frequency pho- 
tons incident on a cloud of electrons which become the seed photons for the 
Comptonization inside the cloud, i.e. the normalization of the distribution 
of multiply (more than once) scattered photons over the number of scatter- 
ings, depends on the plasma temperature. In this section, we shall draw our 
attention to the shape of this distribution, which is also affected, but to a 
lesser degree, by the scattering angular function. 

Sunyacv, Titarchuk (1980) and Payne (1980) obtained analytical formu- 
lae describing the distribution of photons over the time of their escape from 
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a spherical, optically thick plasma cloud. The calculation of these authors, 
carried out in the diffusion approximation, implicitly assumed that the prob- 
abilities are equal for the scatterings in the forward {/ig > 0) and backward 
{/Is < 0) directions. This condition is, for example, indeed met when the scat- 
tering angular function is spherical (p = 1) or Rayleigh {p = 3(1 + /is)/4). 
The latter opportunity is the case when both the photons and the electrons 
are non-relativistic {hu^kT^ <S rrieC?)- However, as we showed in Section 2, 
the forward-backward symmetry is violated if the electrons are substantially 
relativistic. 

It is easy to modify the solution of Sunyaev, Titarchuk (1980) so that 
it would take into account the dependence of the scattering law on plasma 
temperature. To this end we should use the correct value of the diffusion 
coefficient, which is given by (see, e.g., Weinberg, Wigner, 1958) 

^^3(<7/<7t)(1-(/Xs))' ^^^^ 

where cr is the total scattering cross-section, and 

^ _ //isP(/^s) d^j2 

/p(/..)W2 • ^'^^ 

In the non-relativistic limit, a — ut, (/^s) = and D = 1/3, and it is the 
latter value which was used in (Sunyaev, Titarchuk, 1980). 

We shall assume here (as we did in Section 3) that the energy of a photon 
remains small as it diffuses out of the cloud suffering multiple scatterings, 
hence Klein-Nishina corrections are not important. In this case a — (Tt, 
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and (see eq. |T^) only the dependence of (/ig) on the scattering law p{jjs 



determines the magnitude of the diffusion coeffient. 

In the ultra-relativistic case {kTc ^ nieC^), is can readily be found, using 
equation (|^), that (fig) = —1/3 and, consequently, 

f = i (14) 

Therefore, diffusion of low-frequency photons proceeds slower, by a factor of 
4/3, in ultra-relativistic plasma than in cold one. The value (p!^, however, 
has little practical significance. Indeed, a photon acquires a huge amount 
of energy upon scattering from an ultra-relativistic electron of energy 7meC^ 
(z/'/z/ ~ 7^), hence it typically becomes relativistic (being originally low- 
frequency) after few scatterings and our Thomson-limit set-up of the problem 
becomes invalid. 

In the mildly relativistic limit, we derive, using equation (^, that (/is) ~ 
—2kTe/5meC^ and 

^^^^^ " 3{1 + 2kTj5m,c^y ^^^^ 

We would like to point out that our treatment above is similar to that 
in (lUarionov et al., 1979; Grebenev, Sunyaev 1987), where the diffusion 
coefficient that corresponds to the case of scattering of relativistic photons 
by cold electrons was derived. There have also been investigations (Cooper, 
1974; Shestakov et al., 1988) considering the effect of relativistic corrections 
on the so-called transport cross-section, which is used to describe the spatial 
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diffusion of radiation in terms of intensity, rather than in terms of number 
of photons, as we do here. 

Using the diffusion coefficent in the form (|15]), it is easy to recalculate 
distribution functions over the photon escape time that were found by Sun- 
yaev and Titarchuk (1980) for various distributions of the sources of pho- 
tons over a spherical cloud. For example, if the source is situated at the 
center of the cloud, and we are interested in the fate of photons that un- 
dergo more scatterings than the average number, i.e. the dimensionless time 
u = ca-^Nj, ^ Uq = 0.5r^ (where r ^ 1 is the optical depth of the cloud 
along its radius), then 



27r^D(re) 
-T (w) = ^ -— — - exp 

where D{Te) is given by equation (plsl). 



{T + 2D{T,)y 



(16) 



In Fig. 3, we have plotted P{u) for a number of values of kT^. We see 
that the agreement between Monte-Carlo results and the result of equation 



(16) is fairly good up to kTc ~ OArrLeC^. The temperature correction to 
the diffusion coefficient causes photons to spend more time (and suffer more 
inverse- Compton scatterings) in the cloud. 

The knowledge of how photons are distributed over the time they spend in 
a scattering cloud is required in virtually all problems related to Comptoniza- 
tion in astrophysical plasmas. In particular, the shapes of Comptonization 
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spectra and the time delays between the soft and hard components of radia- 
tion coming from variable X-ray sources strongly depend on this quantity To 
solve such problems, it is necessary to consider a related question about the 
evolution of the photon energy with time, which is beyond the scope of this 
letter. We note, however, that much work in this direction (in connection to 
relativistic thermal plasmas) has already been done (Titarchuk, 1994; Hua, 
Titarchuk, 1996; see also references therein). 
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Figure 1: Angular function (in polar coordinates) for the Compton scat- 
tering of photons of energy E by an ensemble of relativistic Maxwellian 
electrons with temperature Tg. Solid lines show Monte-Carlo simulation re- 
sults. Dashed lines represent approximations (when available) given either 
by Eq. @ or by Eq. (H). The dotted lines represent two extreme cases, one of 
vanishing temperature, which corresponds to the Rayleigh angular function, 
and the other of ultra- relativistic electrons {kT^ ^ rrieC^). In the latter case 
(presented only in the left panel), the scattering obeys the law p = 1 — cos a. 
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Figure 2: Albedo of an electron-scattering atmosphere with respect to singly 
scattered photons as a function of the incidence angle a = arccos/xo- The 
radiation is of low frequency, hukTg <^ (meC^)^. Monte-Carlo simulation 
results are shown, (a) The optical thickness of the atmosphere is infinite 

Different curves correspond to various electron temperatures: 



r 



oo 



kTe <^ rrieC^ (solid line), 0.2meC^ (dotted line), O.SmeC^ (short-dashed line), 
2meC^ (long-dashed line) and 3> rrieC^ (dash-dotted line), (b) The same as 
(a), but for a transparent atmosphere with r = 0.5. 
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Figure 3: Distribution of low-frequency photons over the time of their escape 
from a cloud of electrons with an optical depth of r = 5 (i? is the radius 
of the cloud), for various plasma temperatures (curves are labelled with the 
corresponding values of kTe/rrieC^)- The solid lines are Monte-Carlo simu- 
lation results. The dashed lines show the results of the calculation by the 
approximate formula (|TBp, which can be used for describing the exponential 
(power-law in the linear-logarithmic coordinates used) part of the distribu- 
tion function. The diffusion coeffient D appearing in equation (|16D is given by 
Eq. (|T^ in the case kT^ ^ mgC^ and by Eq. (0) in the case kT^ ^ OAnieC^. 
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